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ABSTRACT 

We present analysis of a Chandra observation of the elliptical galaxy NGC 4555. The galaxy 
lies in a very low density environment, either isolated from all galaxies of similar mass or on 
the outskirts of a group. Despite this, NGC 4555 has a large gaseous halo, extending to ~60 
kpc. We find the mean gas temperature to be ^0.95 keV and the Iron abundance to be ^0.5^0. 
We model the surface brightness, temperature and abundance distribution of the halo and use 
these results to estimate parameters such as the entropy and cooling time of the gas, and 
the total gravitational mass of the galaxy. In contrast to recent results showing that moderate 
luminosity ellipticals contain relatively small quantities of dark matter, our results show that 
NGC 4555 has a massive dark halo and large mass-to-light ratio (56.8lg 5 'g MqILbq at 50 
kpc, 42.7^212 at 5r e > lc errors). We discuss this disparity and consider possible mechanisms 
by which galaxies might reduce their dark matter content. 
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1 INTRODUCTION 

The m ajority of galaxies are found in groups and clusters iTulrvl 
Il987l) . and this is particularly true of elliptical galaxies. The 
morphology-density and morphology-radius relations show that el- 
liptical galaxies are most common in the cores of clusters and 
groups lDressleJll980l: iMelnick & Sargendfl977l) . In a hierarchi- 
cal model of structure formation, this can be explained as a prod- 
uctof the processes which form ellipticals. The merger hypothesis 
flbomre & Toomrell 19721) suggests that the product of the merger 
of two spiral galaxies will be an elliptical galaxy. If this is the case, 
then galaxy groups are the most likely location of elliptical forma- 
tion, as these systems have high galaxy densities but relatively low 
velocity dispersions. The merger of groups to form larger galaxy 
clusters naturally leads to a large population of elliptical galaxies 
in the most massive systems. 

Given this model of elliptical galaxy formation, it is per- 
haps unsurprising that the best known elliptical galaxies are found 
in galaxy groups and clusters. Most clusters and many groups 
are dominated by a giant elliptical (or cD) galaxy which lies at 
the bottom of the cluster potential well. X-ray observations have 



shown these objects to be surrounded by haloes of highly lumi 


nous gas l Forman et al. 


1985; 


Trinchieri et al. 


1986). The clus 



ters and groups often have their own haloes of hot X-ray emit- 
ting gas l Kellog g"et alJll975l) . and models of galaxies in the cores 
of such systems suggest that the galaxy halo is probably en- 
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hanced by in flow of ga s from the surrounding intra-cluster medium 
(Mathews & Brishenti 2003). Because of their high X-ray lumi- 
nosities, these galaxies are the most easily observed in X-rays, and 
most detailed analyses of elliptical galaxies focus on th em (e.g. 
iBuote et all2003HSakelliou et al 12001 llones et all2002h . 

However, the location of these galaxies in a dense environ- 
ment, surrounded by a reservoir of high temperature gas means 
that their intrinsic properties must always be in doubt. In particular, 
the question of whether elliptical galaxies produce the observed 
haloes of X-ray emitting gas through stellar mass loss or accre- 
tion from their surroundings is very difficult to answer. It would 
be greatly simplified if galaxies with little or no surrounding intra- 
cluster medium (ICM) could be observed. The importance of this 
issue has recently increased, owing to reports that some ellipticals 
may contain very little dark matter or may lack dark matter haloes 
entirely iRomanowskv et all2003l) . This could provide an explana- 
tion for the long known issue of the large degree of scatter in the 
Lx'.Lb relation for elliptical galaxies - galaxies with dark matter 
haloes have sufficient mass to prevent the escape of hot gas and/or 
accrete more, while those which lack dark matter do not. It is there- 
fore important to observe ellipticals which do not lie at the heart of 
a large group or cluster potential well in order to answer the ques- 
tion of whether ordinary elliptical galaxies can possess significant 
quantities of hot gas, without the aid of aa surrounding deep poten- 
tial well. 

As part of a sample of isolated elliptical galaxies, we have 
used Chandra to observe NGC 4555, a fairly luminous elliptical 
galaxy (log Lb/Lbq = 10J&) at a distance of ~ 90 Mpc. We find 
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R.A. (J2000) 12 35 41.2 

Dec. (J2000) +26 31 23 

Redshift 0.022 kms" 

Distance (H =15) 90.33 Mpc 

1 arcmin = 26.3 kpc 

D25 radius 19.0 kpc 



Table 1. Location and scale of NGC 4555 



that this elliptical, which we show is not the dominant galaxy of any 
group or cluster, has a sizeable X-ray halo. The hot gas in this halo 
can be used to characterise the dark matter halo surrounding the 
galaxy. Throughout the paper we assume Ho=15 km s _1 Mpc -1 
and normalise optical B-band luminosities to the B-band luminos- 
ity of the sun, Lbpi=5.2 x 10 32 erg s -1 . Abundance s are measured 
relative to the ratios of iGreve sse &SairyalNl9 98[). w hich differ 
from the older abundance ratios of Ander s & Grevessd fip89 ) in 
that the solar abundance of Fe is a factor of ~1.4 lower. Details of 
the location and scale of NGC 4555 are given in TableQ 

In Section 121 we give details of the observation and our data 
reduction techniques, and Section[3]contains the results of our anal- 
ysis. Section|4]consists of a discussion of these results, with partic- 
ular reference to the issue of dark matter in early-type galaxies, and 
Section|5|summarizes our results and conclusions. 



N 




20 kpc 



Figure 1. Adaptively smoothed image of NGC 4555 created using 
CSMOOTH with a signal to noise range of 3 to 5. Optical contours are over- 
laid to show the position and extent of the stellar component of the galaxy. 



2 OBSERVATION AND DATA REDUCTION 

NGC 4445 was observed with the ACIS instrument during Chandra 
Cycle 3, Obs ID 2884. A detailed summary of the Chandra mission 
and instrumentation can be found in (Weis skopf et"aHl2002t) . The 
S3 chip was placed at the focus of the telescope in order to take ad- 
vantage of the enhanced sensitivity of the back illuminated CCDs at 
low energies. The instrument operated in faint mode, and observed 
the target for just over 30 ksec. The raw data was reprocessed us- 
ing CIAO v3.0.1 and bad pixels and events with ASCA grades 1, 5 
and 7 were removed. The data were corrected to the appropriate 
gain map, and a correction was made to account for the time de- 
pendence of the gain using the technique described by Vikhlinin 1 . 
A background light curve was produced. Some minor background 
flares were identified and removed, with all periods where the count 
rate deviated from the mean by more than 3a being excluded. The 
effective exposure of the observation after cleaning was 23.3 ksec. 

Background images and spectra were generated using the 
blank sky data described by Markevitch 2 . The data were cleaned 
to match the background, and appropriate responses were created 
using the CIAO tasks MKWARF and MKRMF. As the ACIS instru- 
ments are affected by absorption by material accumulated on the 
optical blocking filter, we applied a correction to the responses. 
When fitting spectra we generally held the absorption fixed at the 
measured galactic value of 1.36x 10 20 cm -2 . Point sources were 
identified using the CIAO WAVDETECT tool with a signal threshold 
of 10—6. We chose to use only the S3 chip for our analysis, as the 
galaxy emission should be entirely contained on S3. This choice of 
signal threshold means that the detection algorithm should identify 
^ 1 false source in the field of view. Once identified, point sources 
were removed from the data, using regions of twice the radius given 
by the detection routine. 

1 http://hea-www.harvard.edu/ alexey/acis/tgain/ 

2 http://asc.harvard.edu/cal/ 



3 RESULTS 

We initially prepared adaptively smoothed images of the galaxy, 
using the CIAO task CSMOOTH. The data were smoothed to show 
features with a signal to noise ratio of 3 to 5. A scaled background 
image was smoothed and removed from the data and correction 
for the variations in exposure across the chip were made. FigureQ 
shows the result of the smoothing with optical contours from the 
Digitized Sky Survey overlaid. From this image, it is clear that the 
X-ray emission extends well beyond the stellar body of the galaxy, 
perhaps as far as 60 kpc in some directions. 



3.1 Two-dimensional surface brightness modeling 

In order to model the X-ray surface brightness distribution of 
NGC 4555 we prepared source and background images in a 0.3-3.0 
keV band, with point sources removed. The energy band was cho- 
sen to focus on soft emission and improve the signal-to-noise of the 
source. Images were binned to a pixel size of 1 ". An appropriate 
exposure map was also generated, and we used the CIAO SHERPA 
package to perform the fitting. As the source image has man y pix- 
els co ntaining few counts (or none), we use the Cash statistic ICashl 
1979) when fitting. This statistic only provides a relative measure 
of the goodness of fit, so that while it allows us to improve fits and 
find the best solution for a particular model, it does not provide an 
absolute measure of the fit quality. We therefore judge whether fits 
are satisfactory (or otherwise) by inspection of azimuthally aver- 
aged radial profiles and residual images. However, the fits are two 
dimensional, and so we can determine parameters such as the ellip- 
ticity of the halo. 

We performed fits using a beta model, de Vaucouleurs model, 
and combinations of the two with a central point source. When 
combining the de Vaucouleurs and beta models, the de Vaucouleurs 
model core radius, axis ratio and position angle were all fixed 
at the optically determined values (r c =16.1l", p.a.=35°, axis ra- 
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Figure 2. Azimuthally averaged surface brightness profile of NGC 4555, 
showing our best fit beta model. The background level is marked by the 
dotted line, the model by a solid line, and the data by error bars. The profile 
was azimuthally averaged using elliptical bins with axis ratio and position 
angle as in the fit. The radius shown is therefore effectively the minor axis 
radius. Note that 1 "=0.437 kpc. 

rcore (") Pfit pos. angle (°) axis ratio 

1.64t°;^ 0.577±0.009 26.77+J-43 1.12±° ° l 

Table 2. Parameters and lcr errors of out best fitting beta model. Position 
angle is measured anti-clockwise from northeast. 

tio=1.26), so that the de Vaucouleurs component would model a 
discrete source contribution distributed in the same way as the stel- 
lar population. We found that the beta model provided an adequate 
fit to the data and inspection of the residual images and azimuthally 
averaged radial profiles did not suggest that more complex models 
improved the fit. We therefore adopt our best beta model fit to de- 
scribe the surface brightness distribution. The fitted parameters and 
given in Table Fi gure|2] shows an azimuthally averaged radial 
profile with the fitted model. 

3.2 Spectral modeling 

As an initial step we chose to fit the integrated spectrum of the 
entire galaxy halo. The spectrum was extracted from a circular re- 
gion centred on the peak of the emission, with radius 150" (~65 
kpc). After removal of regions corresponding to point sources were 
removed, source and background spectra were extracted, appro- 
priate responses created and corrected, and the source spectrum 
grouped to 20 counts per bin. The spectra were fitted using XSPEC 
vl 1.1.0, ignoring energies lower than 0.4 and greater than 8.0 keV. 
The spectrum was fitted fairly successfully with a model consist- 
ing of a 7 keV bremsstrahlung component, and a hot plasma com- 
ponent, modelle d using either the MEKAL iLiedahl et aljfl995l 
iKaastra & Mewdll993h or APEC fcmith et alJl200lh codes. The 
bremsstrahlung component was intended to model emission from 
unresolved point sources. Both bremsstrahlung and power law 
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Figure 3. Integrated spectrum with best fit APEC+APEC+Bremsstrahlung 
model. Lower panel shows residuals to the fit plotted in terms of contribu- 
tions to x 2 ■ 





MEKAL 


APEC 


APEC+APEC 


kT 
kT2 

Zavg 

Si 
Fe 

red. x 2 


0.91±0.04 

awl":" 

0.52±° ;^ 

041 +0.20 

1.093 


0.95±0.04 

0.24«; 49 

0.58t° ;« 

50+ ' 24 
1.056 


0.82±°;° 4 

52+ 3 89 

1 12+ 491 
1 06+ 069 

1.040 


d.o.f. 


162 


162 


160 


Flux 


4.86xl0- 13 


4.75 xlO- 13 


4.54X10- 13 


fbremss 
Lx,gas 


0.186 
4.75 Xl0 41 


0.147 
4.64 xlO 41 


0.126 
4.43 XlO 41 



Table 3. Parameters for our best fits to the integrated spectrum. Temper- 
ature and abundance are given in terms of 90% error bounds. Flux is cal- 
culated for the 0.4-8.0 keV band, is unabsorbed, and is given in units of 
erg s — 1 cm -2 . Gas luminosity is calculated from flux assuming a distance 
of 90.33 Mpc, and is given in units of erg s _1 . The fits were performed 
with A r H fixed at the galactic value (1.36xl0 20 cm -2 ) and included a 7 
keV bremsstrahlung model, fbremss gives the fraction of the total flux orig- 
inating from this component. Z a „ s gives the abundance of all other metals 
apart from Fe and Si. 



models have been successfully used in this role in previous stud- 
ies llrwin et al]|2003l : llrwin et alj|2000h . Results of these fits are 
given in Table|3| 

Although the abundances are rather poorly constrained, free- 
ing Fe and particularly Si greatly improves the fits. A MEKAL 
+ bremsstrahlung model with all metals in solar ratios has a re- 
duced x 2 °f 1-196 for 164 degrees of freedom, a significantly 
poorer fit. Assuming the abundances given by the APEC fit, we 
estimate the masses of Iron and Silicon in the central ~65 kpc 
of the halo to be M Fe =1.49x 1O 5 M and M Si =1.83x 10 5 M© re- 
spectively. This gives an Iron-mass-to-light ratio of only 2.4 x 10 -6 
MQ/LsQ_^considerably lower than the values found for galaxy 
clusters iFinoguenov et all2000l). or even the values suggested for 
small galaxy groups (4.5 X 10~ 4 . lRenzini et alll993h . It seems rea- 
sonable to expect that this relatively small quantity of Iron could be 
produced by the stellar population of the galaxy without external 
enrichment. 

We also attempted to fit the spectrum with more complex 
models, including a two-temperature plasma with bremsstrahlung 
or power law, cooling flow models and multi-temperature models. 
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None of these produced a statistically significant improvement in 
the fit. However, inspection of the spectrum showed that the use of 
a two-temperature plasma model produced a slightly better fit to the 
high energy side of the 1 keV Iron peak, and we therefore record 
this fit in Table|3| Fitting a single temperature model to spectra from 
gas with multiple temperature components around 1 keV is known 
to produce residua ls on both sides of the Iron peak l Bu ote & Fabianl 
ll998tlBuote 2000). so this apparent improvement in fit probably in- 
dicates that NGC 4555 has a multi-temperature halo. It is notable 
that the measured abundances for this fit are considerably higher 
that those found for the single temperature plasma fits, again agree- 
ing with previous studies comparing single- and two-temperature 
fits of multi-temperature gas. However, the abundances are rather 
poorly constrained. To examine whether the large errors in the bins 
above 3 keV might be responsible for this lack of precision, we 
rebinned the spectrum above 3 keV to have at least 20 counts per 
bin after background subtraction. Unfortunately this does not im- 
prove the situation; the errors on temperature and abundance are 
even larger, and the fit has a reduced \ 2 of 1 . 103 (for 85 d.o.f.), not 
a significant improvement on our previous fits. 

We have sufficient detected counts from the galaxy to allow 
us to split the halo into four bins and fit spectra from these indi- 
vidually. We use elliptical annuli with axis ratio and position angle 
taken from the best fitting surface brightness model. The angular 
sizes of the bins are 0-20", 20-45", 45-90" and 90-150", with dis- 
tances measured on the semi-minor axis. The spectra and associ- 
ated responses and background spectra were prepared as described 
above and fitted individually. These spectra were not of the quality 
required to fit individual metal lines, so we used MEKAL mod- 
els with bremsstrahlung components added if necessary. Figure |4| 
shows the fitted temperature and abundance in each bin. In each bin 
we hold the hydrogen column fixed at the galactic value. Only the 
central bin was improved by the addition of a bremsstrahlung com- 
ponent, which contributes ~16 per cent of the emission in this bin. 
A bremsstrahlung contribution is not ruled out in the other bins, 
but provides no significant improvement in the fit. We also note 
that inclusion of such a component does not significantly affect 
the fitted temperature of the MEKAL component. Fit quality in the 
central bin is good, reduced x 2 =0.87 for 29 degrees of freedom. 
Fit quality in the outer bins is considerably poorer, with reduced 
X 2 ~1.3 in bins 2 and 3 and 1.13 in the outermost bin. We tried 
numerous models including MEKAL and APEC plasmas, power 
law and bremsstrahlung components, and multi-temperature com- 
ponents such as CEMEKL and MKCFLOW. None produced any 
significant improvement over a single temperature plasma, and we 
therefore use the parameters determined by these models as the best 
fit. We also tested our background subtraction, using a local back- 
ground extracted from a source-free region of the S3 chip. Bins 1 
and 2 showed no significant change in fit statistic or model param- 
eters when using the local background. The best fit models for bins 
3 and 4 had larger errors on temperature and abundance but were 
consistent with the fits performed using the blank-sky background 
data. Fit statistics for these bins were also slightly altered, presum- 
ably owing to the poorer statistics of the local background data, 
which is the likely cause of the larger error regions. We conclude 
from this that our use of the blank-sky background data is justified 
and gives accurate results. 

Figure|4]shows that temperature rises from the outer parts to- 
ward the core, but turns over and falls in the core. This suggests 
that radiative cooling in the dense core of the galaxy halo is effec- 
tive, and has significantly reduced the mean temperature of the gas. 
The abundance in each bin is more poorly defined, but appears to 
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Figure 4. Best fit projected temperature and abundance for four radial bins 
in the halo of NGC 4555. Grey boxes mark the temperature of each bin, 
and crosses the abundance, with 90% error bounds. The dotted line shows a 
fitted model of the temperature profile, ignoring the central bin, as described 
in Section l33l 

show a general decrease with increasing radius. This is consistent 
with enrichment of the galaxy halo by metals lost from the stellar 
population through stellar winds and supernovae. 



3.3 Mass, Entropy and cooling time 

Given the surface brightness and temperature modelling of the halo 
of NGC 4555, it is possible to estimate three dimensional proper- 
ties such as mass, entropy and cooling time. The density profile 
of the gas can be estimated from the measured profiles and then 
normalised to reproduce the X-ray luminosity of the galaxy, deter- 
mined from our best fitting APEC model. Given the density profile 
we can use the well known equation for hydrostatic equilibrium, 

kTr / (Mpgas ^ dlnT^ 
firripG \ dlnr 



M to t{< r) = 



(1) 



d\ur ) ' 

to calculate the total mass within a given radius. From as density 
and total mass, we can calculate parameters such as gas fraction, 
cooling time and entropy where entropy is defined to be 

T 



S 



(2) 



We estimate the errors on the derived values using a monte- 
carlo technique. The known errors on the temperature and surface 
brightness models, and on other factors such as the total luminosity, 
are used to randomly vary the input parameters. We then generate 
10000 realisations of the derived parameters profiles, and use these 
to calculate the la error on each parameter at any given radius. 

The only issue which arises in these calculations is the ques- 
tion of how well we can model the temperature profile, which de- 
clines sharply in the galaxy core. The central bin has no effect on 
the value of the parameters at larger radii and increases the com- 
plexity of the required model, and we therefore choose to exclude 
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the central temperature bin, and model the temperature profile as 
if there were no central cooling. The remaining three bins can be 
well described with a quadratic, which is shown as a dotted line in 
Figure|4] Based on this model of the temperature we calculate the 
parameters shown in Figure|5| 

In order to calculate the mass-to-light r atio, we use the H band 
near infra-red surface brightness profile of iGavazzi et alj |2000). 
This profile is made up of a bulge component described by a de 
Vaucouleurs profile, with r e =4.0l", and an exponential disk with 
r e =17.42". We assume that these measurements hold for the B 
band optical light, and normalise the profile to the B band lumi- 
nosity, with a bulge-to-disk ratio of 0.33. We note that the mass-to- 
light ratio at the inner limit of the plot is ~9, a little higher than 
the value of 5-8 generally assumed for the mass-to-light ratio of the 
stars alone. 

The outer limit of the plots is determined by the radius to 
which we can measure the temperature, and the inner limit by the 
radius of the innermost temperature bin. The inner limit is chosen 
to exclude the 5 kpc radius region in which we know our tempera- 
ture model to be wrong. We also note that both the temperature and 
surface brightness model profiles were determined using the semi- 
minor axis as the radius descriptor. The profiles shown in Figure|5| 
should therefore be considered to show properties which are az- 
imuthally averaged around ellipses whose position and ellipticity 
are determined by the best fitting surface brightness model, with 
the radial scale indicating the semi-minor axis. 



3.4 Point sources 

As described in Section|2| we used the CIAO WAVDETECT tool to 
identify point sources in the field of view. Several sources are found 
to lie within the extended emission surrounding NGC 4555, but 
only two lie within optical extent of the galaxy, defined by the D25 
ellipse. Of these, one is coincident with the galaxy core. Surface 
brightness fitting of the galaxy X-ray halo (see Section l3~Tl does not 
show strong evidence for a central point source, and we conclude 
that the identification from WAVDETECT actually corresponds only 
to the peak of the galaxy halo, not to a separate source. 

The remaining source lies ~15" NE of the galaxy centre, and 
does not correspond to any feature visible in the DSS optical or 
2MASS infra-red images. There are no objects listed at its posi- 
tion in NED. We extracted a count rate for the source, in the full 
Chandra band and in a 0.5-2.0 keV band. Background subtraction 
was carried out using a region immediately surrounding the source, 
and should therefore account for contamination by the galaxy halo. 
Based on the background subtracted count rate we estimated the 
flux and luminosity of the source assuming a bremsstrahlung model 
with kT=5 keV, and a power law model with T= 1 .96, typical of high 
luminosity point sources in other galaxies llrwin et alj|2003l) . The 
results are given in Table|4| 

If this source is in fact part of the NGC 4555 system and is ac- 
tually a single object rather than an unresolved cluster of sources, 
then it is extremely luminous, log Lx~39.43 ergs -1 . We can 
estimate the probability of finding a background source with the 
measured flux within the D25 ellipse based on the sources found 
in the Chandra deep field south flbzzi et al]|200lh . The D25 el- 
lipse has an area of ~3.62xl0~ 4 square degrees, which means 
that we would expect to find ~0.02 sources with the observed (or 
greater) flux in that area. Assuming that the background distribu- 
tion of sources is the same as the Chandra deep field, this suggests 
that the source is probably part of NGC 4555, and must therefore 
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D A 
K./\. 


7/1 

1Z. J J .4 1 . /4 


R.A. error (") 


±0.127 


Dec. 


+26:31:35.2 


Dec. error (") 


±0.123 


Counts (0.1-12.0 keV) 


15.04 


Counts (0.5-2.0 keV) 


15.76 


Flux (BR, erg s _1 cm -2 ) 


2.719X10" 15 


Flux (PL, erg s" 1 cm -2 ) 


2.768x10" 15 



Table 4. Properties of the off centre point source in NGC 4555. Fluxes are 
calculated assuming no absorption, a 0.5-2.0 keV energy band and either a 
5 keV bremsstrahlung (BR) model or a power law (PL) model with T=l .96. 

be classed either as an ultra-luminous X-ray source (ULX), or con- 
sidered to be an unresolved cluster of sources. 



4 DISCUSSION 
4.1 Environment 

NGC 4555 was selected as part of a sample of isolated ellipti- 
cals, extracted from the Lyon-Meudon Extragalactic Data Archive 3 
(LEDA). Sample galaxies were selected using the following crite- 
ria: 

(i) Morphological type T ^ —3, i.e. early-type galaxies. 

(ii) Virgocentric flow corrected velocity v ^ 9000 km s _1 . 

(iii) Apparent B-band magnitude Bt 14.0. 

(iv) Galax y not listed as a member of a Lyon Galaxy Group 
(LGG. lGarciall993h . 

The restrictions on apparent magnitude and recession veloc- 
ity were imposed to minimise the effect of incompleteness in the 
catalogue. The LEDA catalogue is known to be 90 per cent com- 
plete at Bt = 14.5 lAmendola et alll 997^. so our sample should 
be close to 100 per cent statistically complete. The selection pro- 
cess produced 330 galaxies which could be considered as potential 
candidates. These were compared to the rest of the catalogue and 
accepted as being isolated if they had no neighbours which were: 

(i) within 700 km s _1 in recession velocity, 

(ii) within 0.67 Mpc in the plane of the sky, 

(iii) less than 2 magnitudes fainter in Bt- 

The criteria were imposed to ensure that the galaxies did not 
lie in groups or clusters, and to ensure that any neighbouring galax- 
ies were too small to have had any significant effect on their evolu- 
tion or properties. 

To check the results of this process, all candidate galaxies 
were compared to the NASA Extragalactic Database (NED) and 
the Digitised Sky Survey (DSS). A NED search of the area within 
0.67 Mpc of each candidate identifies galaxies not listed in LEDA. 
We also examine DSS images of this region for galaxies of similar 
brightness to the target which are not listed in either catalogue. The 
process produced 40 candidate isolated elliptical galaxies, of which 
NGC 4555 is one. 

Although NGC 4555 meets the criteria described above, it 
does have a number of galaxies relatively close to it. We have 
therefore examined the surrounding galaxy population in order to 
determine whether NGC 4555 is on the outskirts of a group or 
cluster. Figure [6] shows a DSS image of the region surrounding 

3 http://leda.univ-lyonl.fr 
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Figure 5. Projected temperature and deprojected gas density, gravitational mass, gas mass, gas fraction, entropy, cooling time and mass-to-light ratio for 
NGC 4555. The inner boundary of the plots is 5 kpc, the radius within which we know our temperature model to be inaccurate. Solid lines show values derived 
from the best fitting temperature and surface brightness models, grey regions show lcr errors. 



NGC 4555, which is marked in the centre of the image. Scale is 
shown by the two large circles, which have radii of 0.67 and 1 
Mpc. All galaxies listed in NED within 1 Mpc which have mea- 
sured recession velocities are marked on the plot. Crosses indicate 
galaxies >700 km s _1 away from NGC 4555, circles those within 
this velocity range. All galaxies marked by circles have apparent 
magnitudes at least 2 magnitudes fainter than NGC 4555, with the 
exception of IC 3585, a nearby SO galaxy which NED shows to be 
-600 kpc and 687 km s" 1 from NGC 4555. It is only ~1.5 mag- 
nitudes fainter than our target, and so would appear to violate our 
isolation criteria. However, the velocities available from NED have 
not necessarily been corrected in the same way, so the LEDA veloc- 
ities should provide the more accurate measure of relative distance. 
Using the LEDA Virgocentric flow corrected velocities, the differ- 
ence between the two galaxies is 758 km s~ , putting IC 3585 just 
outside our chosen velocity limit. 

Figure Q shows histograms of the local velocity field within 
1 Mpc of NGC 4555. In the upper panel we show galaxies found 
in NED, which were required only to have a measured redshift. 
In the lower panel we show galaxies found in LEDA, which were 
required to have a measured redshift and apparent magnitude. We 
have marked galaxies whose magnitudes indicate that they have a 



B-band luminosity in excess of 10 10 erg s _1 .The remaining galax- 
ies have 2 x 10 9 < L B < 10 10 erg s _1 . NGC 4555 is represented 
in each plot by the shaded region at lowest recession velocity. 

Although the plots do suggest some sort of extended structure, 
it is clear that the galaxies in this region do not form a relaxed 
group. Similarly, it is clear that NGC 4555 is separated from the 
other galaxies of significant mass (for which Lb is a proxy) by 
significant distance and velocity difference. It seems possible that 
NGC 4555 is part of a small filament of galaxies, or that it (with 
some small companions) is falling into a very poor group to the 
north. It is very unlikely that NGC 4555 is the dominant galaxy of 
a virialised group, and it therefore must be assumed that it does not 
lie in the centre of a massive group-scale dark matter halo. 

NGC 4555 does not appea r in the sample of "very isolated 
early-type galaxi es" defined by Stocke et alj 120041) . This sample 
is drawn from the Karache ntsev; " ^973lU solated galaxy catalogue, 
which was compiled by examining images of the environment of 
all bright (m_g ^ 15.7) galaxies in the Zw ickv et alj 1 1968) cata- 
logue. Galaxies were classified as isolated if they had no compan- 
ion within 20 galaxy diameters, where a companion was defined as 
any galaxy with angular diameter within a factor of 4 of that of the 
candidate. Stocke et al have further cleaned the resulting sample by 
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Figure 6. A Digitized Sky Survey image of the environment of NGC 4555 
with positions of galaxies overlaid. The two large circles mark regions of 
radius 0.67 Mpc and 1 Mpc, centred on NGC 4555. Small circles mark 
galaxies whose recession velocities differ from that of NGC 4555 by <700 
hns -1 , crosses those >700 kms - 1 away. For further details, see the 
text. North is to the top in this image and west is to the right. The linear 
feature running across the lower third of the figure is the junction between 
DSS plates. 



examining the Palomar Observatory Sky Survey (POSS) plates to 
confirm isolation, and performing deeper imaging around a num- 
ber of the galaxies. This should lead to a sample with few galaxies 
falsely identified as isolated, but will rule out NGC 4555 because 
of similarly sized galaxies nearby which are in fact distinct from it 
in redshift space. It is possible that NGC 4555 is ruled out by the 
presence of NGC 4565, the large edge on spiral galaxy at the lower 
edge of Figure|6| which has a recession velocity >5000 km s _1 
smaller than NGC 4555. 

We also note that NGC 4555 has been identified as a member 
of a galaxy group by two previous studies. iMahtessiar] ll998l) sug- 
gest that it is part of a triplet with IC 3582 (which is fainter than our 
magnitude cut off) and NGC 4556 (whic h has a recession ve locity 
~800 km s" 1 than that of NGC 4555). iRamella et alj Jiool use 
a friends of friends algorithm to search for associations in the Up- 
dated Zwicky Catalog (Falco et aljfl999h and Southern Sky Red- 
shift Survey Ida Costa et al. 1998) galaxy catalogues. They suggest 
that NGC 4555 is part of a group of 7 objects with a range of re- 
cession velocities (6615-7819 kms -1 ). They calculate the group 
velocity dispersion to be 568 km s" 1 and from this estimate the 
virial radius (1.64 Mpc), total mass (3.7xl0 14 Mq) and mass-to- 
light ratio (~1740 MqILbq)- These values seem rather large for 
a system so poor in galaxies, and are particularly surprising given 
the lack of a group X-r ay halo. Fro m the mass -tempe rature relation 
of clusters and groups ISanderson & Ponmanll2003T) we would es- 
timate that a system of this mass would have an X-ray temperature 
of kT~8 keV, and it seems unlikely that such a large dark matter 
potential could have failed to accrete some gas during its formation. 
We therefore consider it likely that the identification of NGC 4555 
and its neighbours as a collapsed group is faulty, and that they at 
best form a loose association, not dominated by a single large dark 
matter halo. 
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Figure 7. Redshift distribution of the galaxies within 1 Mpc (in the plane 
of the sky) of NGC 4555. The upper panel shows galaxies with redshifts 
listed in NED, while the lower panel shows galaxies listed in LEDA with 
both redshifts and B-band magnitudes. Shaded areas show galaxies whose 



B-band luminosity is > 10 ergs 



4.2 Comparison with other systems 

As mentioned in Section^ recent studies of the velocity dispersion 
profiles of some "ordinary" ellipticals have s hown them to have 
lower mass-to-light ratios than were expected jRomanowskv et all 
2003). This suggests that the dark matter content of these galaxies 
is minimal, at least within the radius to which the measurements 
extend (~5r e ). The profiles are incompatible with the dark matter 
profiles predicted by simulations of structure formation in a cold 
dark matter dominated universe, leading to the conclusion that ei- 
ther the haloes take the predicted form but have much lower masses 
than expected (the results are consistent with zero dark mass), or 
that the dark matter halo takes a very different form, perhaps with 
most of the mass at larger radii. Four elliptical galaxies, NGC 821, 
NGC 3379, NGC 4494 and NGC 4697, have been shown to have 
this unexpected lack of dark matter. Of these, three are in small 
groups and are the brightest elliptical in each system. NGC 821 
is a rather more isolated system, not part of any known group or 
cluster. All of the galaxies are relatively X-ray faint, undetected (or 
only marginally detected, in the case of NGC 4697) by ROSAT, and 
two have bee n shown to possess relatively sm all amounts of X-ray 
emitting gas llrwin et alKOOCtlO' Sullivan & Poninar]|2004 . This 
lack of a sizeable X-ray halo could be a consequence of the lack 
of dark matter in these galaxies, as the stellar mass alone would be 
insufficient to retain a large gaseous halo. 

For comparison with NGC 4555, we use the mass-to-light 
rat ios at 5r e for NGC 4494 NGC 3379 and NGC 821, quoted 
by iRomanowskv et all 120031) . An overall mass-to- light ratio for 
NGC 4697 of M/L=11M Q /L S was found by iMendez et ail 
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Galaxy Radius M/L 

n (m q /l bq ) 

NGC 821 5r e =250 13-17 
NGC 3379 5r e =175 5-8 

NGC 4494 5r e =245 5-7 

NGC 4555 5r e =80.6 43~6 

5r eb =20.1 13.8 

5r ed =87.1 46.4 

Table 5. Mass-to-light ratios for NGC 4555 and three ellipticals believed 
to have little or no dark matter. r e is the effective radius of the galaxy as a 
whole, r e t and r e( j are the effective radii of the bulge and disk components 
of NGC 4555. 



1200 ill , but this is for a radius of ~3r e . For NGC 4555 itself we 
have three different measurements of r e , one overall value, which 
assumes a de Vaucouleurs profile, and th e effec tive radii of bulge 
and disk components from lGavazzi et all l2000h . who decompose 
the surface brightness profile of the galaxy into two components. 
Their result is rather confusing, as their profiles assign only 33 
per cent of the H band light of the galaxy to the bulge compo- 
nent, which would suggest that NGC 4555 is a misclassified SO. 
However, NGC 4555 is classified as an elliptical in both NED and 
LEDA, and in fact Gavazzi et al. list its morphological type as ellip- 
tical. We therefore believe that the overall effective radius probably 
gives the best indication of the scale of the galaxy, but quote mass- 
to-light ratios at three radii. The ratios are listed in Table [5] Using 
the overall effective radius, it is clear that NGC 4555 has a consid- 
erably larger mass-to-light ratio than the three other ellipticals. The 
Gavazzi et al disk component effective radius produces a similar re- 
sult, and only if the bulge component is used do we find comparable 
figures. We assume that this is a poor measure of the true scale of 
the NGC 4555, and that the galaxy has a larger mass-to-light ratio 
than the three Romanowsky et al. ellipticals. late 

Another class of elliptical galaxies which should be compared 
with NGC 4555 are the fossil groups. These are systems in which 
all the major galaxies of a group collapse at an early epoch, merg- 
ing to form a single giant elliptical embedded in a group scale dark 
matter and X-ray halo l Ponman & Bertradfl993l) . Because these 
systems appear as a single elliptical with no neighbouring galaxies 
of similar size, our optical selection criteria would likely identify 
them as candidate isolated ellipticals. Fossil groups can be identi- 
fied by four main features IJones et alj|2003h . These are: 1) They 
have an X-ray luminosity typical for a galaxy group (Lx > 10 42 
erg s _1 ), 2) Their X-ray halo is highly extended, 3) The domi- 
nant elliptical is at least 2 magnitudes brighter than the other group 
galaxies, and 4) The dominant elliptical is surrounded by a halo 
of faint galaxies, the unmerged members of the galaxy group. Cri- 
teria 3 and 4 can be refined in that a luminosity function of the 
galaxies in the fossil group will have only one galaxy with opti- 
cal luminosity greater than L*, the luminosity of this galaxy will 
be significantly higher than would be expected from the luminos- 
ity functions of other groups, but the tail of dwarf galaxies will be 
relatively similar t o that of other galaxy groups (see in particular 
Fig. 7, IJones et aljkood) . NGC 4555 fails the first of these cri- 
teria in that its X-ray luminosity is only ~4.4-4.7x 10 41 erg s _1 
(model dependent). The extension of the halo of NGC 4555 is also 
smaller than any known fossil group. The one possible exception is 
NGC 6482, a fossil group observed with Chandra, whose halo ex- 
tends off the ACIS-S detec tor and so is as yet poorly characterised 
as regards extent (Khosroshahi et al. 2004). The optical luminosity 
and isolation of NGC 4555 has been dealt with in Section l4~T1 and 



from this the galaxy appears to be isolated enough to meet condi- 
tion 3 for fossil group status. Testing condition 4 is more difficult, 
as we do not have redshifts for the faint galaxies around NGC 4555. 
However, we extracted the positions of all galaxies without redshift 
listed in NED within 0.444° of NGC 4555 (equivalent to 0.7 Mpc 
at the distance assumed, 90.33 Mpc) and plotted a radial profile of 
surface number density of galaxies around our isolated elliptical. If 
NGC 4555 were a fossil group we might expect to see higher num- 
ber densities around it, caused by the surrounding halo of fainter 
group members. We note that the galaxies listed have magnitudes 
as faint as 19.5 mag. and we are therefore likely to be missing some 
of the faintest galaxies which have not been identified and listed 
in NED. However, we find no evidence of an overdensity of faint 
galaxies around NGC 4555 and it seems unlikely that the inclusion 
of a small number of fainter objects could change this result. This 
combination of environmental and X-ray properties argues strongly 
against NGC 4555 being a fossil group. 

It is also possible to compare our results with other X-ray 
mass estimates from the literature. A number of early-type galax- 
ies have mass estimates available, but the requirement for high 
quality data means that the majority of these galaxies are highly 
luminous and reside in the cores of groups and clusters. Tech- 
niques similar to ours have bee n used to produce mass profiles 
of galaxies suc h as NGC 507 <Paolillo et alj|2003h . NGC 1399 
and NGC 1 404 iPaolillo et alJl2002h . NGC 2563. NGC 4325 and 
NGC 2300 jMushotzkv et al]bo03h. NG C 4472, NGC 4636 and 
NGC 5044 (Mathews & Brishenti 2003, and references therein). 
However, these are all group or cluster dominant galaxies, with 
the exception of NGC 4472, which dominates a subclump in the 
Virgo cluster, and NGC 1404, which has a truncated halo, probably 
caused by intera ction with Fornax clus ter gas. Using an alternate 
technique. lLoewenstein & W hite ( 1999) determined constraints on 
mass and mass-to-light ratio for a sample of ~30 galaxies for which 
global temperature measurements are available, but more detailed 
density and temperature profiles are not. The authors find that the 
relation between X-ray temperature and optic al velocity d isper- 
sion determined from the sample of iDavis & White! Il99rjl im- 
plies a fairly constant mass-to-light ratio within 6r e of ~23 /175 
Mq/Lbq- Once again, however, the sample of galaxies used is 
heavily weighted toward those ellipticals which dominate groups 
and clusters - 20 of the 30 galaxies listed by Davis & White are 
group or cluster dominant ellipticals, and a further 3 should prob- 
ably be excluded - M32, a dwarf elliptical interacting with M31, 
NGC 4406, undergoing strong interaction with the Virgo ICM, and 
NGC 4472 which is mentioned above. It is also worth noting that 
of the remaining 7 galaxies, 6 are located in either the Virgo or For- 
nax clusters, and could have had their dark matter halos altered by 
interactions with their environment. 

A third technique for determining the presence of dark mat- 
ter from X-ray observations has been demonstrated i n a series of 
papers studying NGC 1332, NGC 3923 and NGC 72 jBuote et alJ 
l2002t iBuote & Canizareslll998t Il997l Il996l Il994l) . These galax- 
ies have non-spherical stellar and X-ray distributions, and the au- 
thors employ geometric arguments to show the need for a dark 
matter component to explain their ellipticity and the relative po- 
sition angles of the gas and stars. They are able to model the dark 
halo and produce strong constraints on the total mass and mass- 
to-light ratio, under the assumptions of hydrostatic equilibrium 
and that rotation has a minimal influence on the potential. These 
models suggest that NGC 720 has a mass-to-light ratio of ~19 
Mq/Lbq at 3r e , NGC 3923 has M/L=17-32 MqILbq at ~47 
hrs kpc, and NGC 1332 has M/L-31-143 MqILbq within ~46 
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hrs kpc. However, these galaxies are all group dominant ellipticals, 
and while the groups are not as X-ray luminous as objects such as 
NGC 5044, they are gravitationally bound systems, as confirmed 
by the group-scale X-r ay haloes around NGC 720 and NGC 3923 
iMulchaev et all2 003). The mass-to-light ratios for these galaxies 
are in general agreement with those found by the other methods de- 
scribed above, suggesting that at the radii of interest to us (~5-6r e ), 
most ellipticals for which X-ray estimates are available have mass- 
to-light ratios of at least 20 Mq/Lbq, and in some cases consider- 
ably higher. At larger radii, the mass-to-light ratio can rise consid- 
erably higher, above 100Mq/Lbq, indicating the influence of the 
group or cluster dark matter halo. 

The results of these X-ray mass studies show that dark matter 
is present in and around the galaxies observed. This is an important 
result, but the fact that these galaxies are almost all in the cores 
of much larger systems must raise the issue of whether the dark 
matter profiles derived for them are really describing the galaxy, the 
group halo, or some combination of the two. The galaxies for which 
individual mass profiles or determinations are available can tell us a 
great deal about ellipticals at the centres of larger structures, but we 
cannot know what the influence of their environm ent is. Ther e are a 
small number of objects in the Loewenstein & White 1 1999) study 
which might be individually useful, but unfortunately the technique 
used applies to the full sample of galaxies and does not provide 
mass-to-light ratios for each elliptical. We therefore conclude that 
while we can compare our results for NGC 4555 to these X-ray 
mass estimates, we must do so cautiously, considering that we may 
be comparing systems of quite different scale and content. 

Two conclusions might be drawn from these comparisons. 
Firstly, as we have already demonstrated, NGC 4555 is not in the 
core of a virialised group, and is therefore not surrounded by a 
group or cluster scale dark matter potential. Its observed proper- 
ties therefore confirm that elliptical galaxies can possess dark mat- 
ter haloes of their own, regardless of their environment. The mass 
of dark matter seems to be comparable to that found for ellipti- 
cals in he cores of groups and clusters. It is to be expected that at 
larger radii, group and cluster dominant ellipticals would have con- 
siderably higher mass-to-light ratios that NGC 4555, and it would 
be interesting to extend our mass profile further to investigate this. 
The low dark matter content found in the three ellipticals studied 
by Romanowsky et al. show that while ellipticals can possess dark 
haloes, not all of them do. This raises an important question; what 
determines whether elliptical galaxies have dark matter haloes? 

One possibility is that all elliptical galaxies are formed with 
dark matter haloes, but some later lose them through interactions 
with other galaxies. Close interactions between galaxies can cause 
tidal stripping of the dark matter halo I Mat hews & Brighentill997() 
as well as gas and stars. Simulations of interactions among multiple 
galaxies in a compact group suggest that a large amount of the dark 
matter may in fact be dispersed, forming a common halo but not 
bound to any particular galaxy lBarnesll989l) . The simulations also 
show that the dominant galaxy of the group retains a sizeable dark 
matter halo, but it is possible that there are circumstances in which 
this would not be the case. NGC 4494, NGC 3379 and NGC 4697, 
all of which are the most luminous and presumably most massive 
elliptical in their groups, might have lost their dark matter in this 
way. 

NGC 821 is more difficult to explain as a product of tidal strip- 
ping. The galaxy is relatively isolated, not a member of any group 
or cluster, and in fact meets the isolation criteria described in Sec- 
tion l4.ll As it has no massive neighbours and does not appear to be 
part of a larger structure, there is little chance that it has suffered 



tidal stripping. Another method of removing its dark matter halo, 
or forming such a galaxy without a halo, therefore appears to be 
required. 

Whereas we can compare the mass and mass-to-light ratio 
of NGC 4555 with the galaxies studied by Romanowsky et al. , 
comparison of the gas mass and gas fraction with other systems is 
hampered by the lack of detailed studies of bright ellipticals out- 
side the cores of groups and clusters. As mentioned previously, al- 
most all elliptical galaxies for which mass and gas mass profiles 
have been calculated are the dominant galaxies of larger structures. 
Their properties are likely to be affected by the surrounding po- 
tential, and models of galaxies embedded in a dense intra-group 
medium show that the galaxy can accrete gas from its environ- 
ment iBrighenti & Mathews 1999), thereby changing its gas frac- 
tion. An alternative is to compare NGC 4555 to a sample of poor 
groups. These have gas temperatures similar to that we observe in 
NGC 4555, and sufficient numbers have been studied to make sam- 
ples re liable. We use the 0.3-1.3 keV sample of ISanderson et all 
(2003) which contains two elliptical galaxies but is dominated by 
poor groups. 

For an accurate comparison, it is necessary to compare prop- 
erties at a common radius, relative to the overall size of the system. 
This is usually done using R200 (where -R200 is a good approxi- 
mation of the virial radius). Unfortunately we do not have a mea- 
sured value of -R200 for NGC 4555 and we cannot calculate one 
based on our three dimensional models, as the temperature model 
is not accurate beyond the outer radius we have used. Based on the 
Sanderson et al. groups, a typical -R200 for a system of this tem- 
perature might be ~500 kpc. If we adopt this value, we see that 
while the 0.3-1.3 keV poor groups have gas fractions of ~1 per 
cent at 0. lxifeoo, NGC 4555 has a gas fraction a factor of ~5 
lower. This result is of course dependent on the value of -R200 cho- 
sen. The Sanderson et al. systems have values of R200 ranging from 
~200-800 kpc, but even if we assume a larger R200 ( and therefore 
measure gas fraction at a larger radius) we still find that NGC 4555 
has a gas fraction considerably lower than that of poor groups. 



5 SUMMARY AND CONCLUSIONS 

We have used Chandra to observe the relatively isolated early-type 
galaxy NGC 4555. An examination of its environment suggests that 
the galaxy is not a member of a virialised group, though it may be 
part of a loose association or filament of galaxies. It is unlikely 
that the galaxy is in the core of a larger, group-scale dark matter 
potential, as is the case with many of the early-type galaxies whose 
X-ray properties have been studied in more detail. Despite the lack 
of a surrounding group potential, we find that the galaxy possesses 
an extended gaseous halo with a temperature of kT~0.95 keV and 
Fe abundance ~0.5Zq. 

We measure the surface brightness distribution of the gaseous 
halo and find that it is reasonably well described by a single beta 
model, supporting the spectral results which suggest that emission 
from gas dominates over emission from point sources. We also 
measure the temperature and abundance profiles and find evidence 
for a central cooling region, though confirmation of this would re- 
quire fitting a multi-temperature model which at present we have 
insufficient counts to do. Assuming the gaseous halo to be in hy- 
drostatic equilibrium, we use these measurements to estimate the 
gas mass, entropy, cooling time, total gravitating mass and mass-to- 
light ratio of the system. We find a mass-to-light ratio of 42.7^J,i '% 
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Mq/Lbq at 5r e , demonstrating that dark matter makes up an im- 
portant part of the galaxy mass budget. 

A recent optical study of the dark matter content of three ellip- 
ticals by Roman owskv et all 120031) shows a quite different result. 
All three galaxies have very low dark matter content, and two are 
consistent with having no dark matter at all, at least out to 5r e . At 
least one of these galaxies (NGC 4494) is X-ray faint, and two of 
them are members of galaxy groups. We suggest that the X-ray lu- 
minosity of early-type galaxies, which is dependent on the size of 
X-ray halo which they can maintain, may be an indicator of their 
possession, or lack of, a dark matter halo. This raises the question 
of how galaxies could lose (or gain) an extensive dark matter halo, 
and how we might distinguish between different processes which 
could affect such haloes. This is clearly a question which deserves 
further consideration, and it seems likely that both improved mod- 
els and further observations of early-type galaxies will be required 
before it can be answered. 
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The isolated elliptical NGC 4555 observed with Chandra 
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